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ABSTRACT 

The radiated acoustic field from a point-driven, 
fluid-loaded plate was determined several years ago. 
However, current operational problems’ force a re- 
examination in the light of certain acoustic anomalies. 

The types of waves that propagate in a solid plate 
are reviewed for their relative importance. 
Experimental models were constructed in an anechoic 
acoustic tank in an attempt to model the aforementioned 
anomalies. 

It is found that the effect of a plexiglass plate 
placed at the boundary of a method-of-images geometry is 
that of a phase delay only. Therefore, mode conversion 
from compressional waves to flexural waves is thought to 
be insignificant. 


A second experiment is described in which a plexi- 


glass plate was driven with a small longitudinal 
vibrator. Experimental radiation patterns suggest that 
of a quasi-piston. Plexiglass is an acoustically slow 


media, therefore the coincidence effect is not possible. 
Further, its absorption of acoustic waves make it a poor 
choice’ for experiments where wave propagation is 


important. 
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10 


Pm 
OPs=t Pp, 
o 


r 


6,6(x,yY,X,t) 


> + 
we W(X,Y,Z,t) 


Wy Wa 


LIST OF SYMBOLS 
(GREEK ALPHABET) 


Fluid loading factor 
Plate shear parameter 
Loss factor 

Angle, defined locally 
Wavelength, wavelength in fluid 
Wavelength at coincidence 
Flexural wave length 

Density of fluid 

Density of plate 

Density per unit area 

Poisson's ratio 


Phase delay 


Scalar displacement potential, azimuthal 
angle 


Vector displacement potential 


Radian frequency, radian coincidence 


frequency 


Vector Operations 


The gradient 


The LaPlacian 
The curl 


The Biharmonic operator 


Lak 


I. INTRODUCTION 


Ae THESIS INTROEUCI (TC) 

The experiments on which this thesis is based were 
attempts to physically model a sound source mounted on 
the ice in the Arctic using a sheet of plexiglass and 
piezoelectric transducers in an anechoic water tank. 
Specifically, the sound field radiated into the fluid 
was of interest. 

There are many wave types that propagate in a solid 
and along a solid-fluid interface. Each was examined 
for its relative importance, so that attention could be 
concentrated on those that are more important to the 
radiation of sound into the water from a sound source 
mounted on the ice. The Ewing and Press’ [Ref. 1] method 
of solution source wave propatation in ice far from the 
source. Stein [Ref. 2] considered radiation from sound 
sources mounted on the ice. 

Two experiments were performed. The first was an 
investigation on the effect of a plate placed at the 
free boundary of the classic method-of-images (Lloyd's 
mirror) geometry. Waves emanating from a source near a 
pressure-release boundary will develop an interference 
pattern with waves reflected from the boundary that can 


be detected by a receiver. It is hypothesized that a 
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plate placed at the boundary will affect the inter- 
ference pattern. 

In the second experiment performed, the source was 
attached to the ice in an attempt to duplicate the 
“acoustic mode," named by Stein [Ref. 2, p.100], of the 
Arctic ice sound propagation. The theoretical develop- 
ment of the Chapter uses the Mindlin-Timoshenko thick 
plate equations to explain the effects of coincidence 
and fluid loading on the radiated sound field froma 


point-driven, fluid-loaded plate. 


B. BACKGROUND 

The classical two-dimensional theory of flexural 
waves in elastic bars in vacuo, using Lagrange's 
equation, was solved in 1889 by Lord Rayleigh under the 
assumption that the wavelengths were long with respect 
to the thickness of the plate. Thus, the equations are 
invalid for high frequency waves or for sharp 
transients. There are two effects that were neglected: 
rotatory inertia and transverse shear deformation. 
Later, however, Rayleigh introduced into the equations 
the effect of rotatory inertia, but because transverse 
shear deformation was not included, this correction was 
too large. Timoshenko [Ref. 3] included the effect of 
transverse shear deformation ina study of flexural 


waves of Dars. 


Mindlin [Ref. 4] in a study of flexural waves in 
plates included the effects of both rotatory inertia and 
shear deformation in his now classic paper in which are 
developed the so-called Mindlin-Timoshenko thick plate 
equations. Dyer [Ref. 5], using the Mindlin-Timoshenko 
equations, showed the effects of transverse shear to be 
important even for thin plates. Skudrzyk et al [Ref. 6] 
developed a theory of flexural waves ina driven plate 
Similar to that of the theory of reverberation in room 
acoustics. 

In 1966, David Feit [Ref. 7] calculated the acoustic 
pressure radiated into a fluid by an infinite plate 
driven by a point force, using the Mindlin-Timoshenko 
thick plate equations. His solution used the’ saddle 
point method of steepest descent to evaluate the aametn 
sure integral. Stuart [Ref. 8] explored the solution of 
the Mindlin-Timoshenko equations with then-new insights 
concerning the so-called leaky waves emanating from the 
plate. He asserted that previous discrepancies in the 
saddle point technique solution (a coincidence frequency 
lower than predicted, and small but unaccounted for 
lobes in the far-field pattern below coincidence) vaias 
due to neglected poles occurring near the saddle point, 
called leaky wave poles, and that these influence the 


far-field pressure. 


Creighton [Ref. 9] contradicted these findings, 


demonstrating that ". . . the far-field contributions 
associated with .. . the leaky waves may be altered 
merely by a change in the integration path. . ." and 


that the contributions decayed exponentially with depth 
into the water and did not represent physically meaning- 
ful free modes, except for subsonic surface waves. 

As the theory of waves in a plate developed, so did 
the theory of wave propagation in the ice. The theory 
of propagation of elastic waves ina floating ice sheet 
was given by Press and Ewing [Ref. 1], by considering 
the ice as an infinite plate and using classical wave 
theory. Press and Ewing [Ref. 10] then extended their 
results to include the effects of the air. A further 
discussion of their work as a whole is Pound er ath: , 
Jardetzky, and Press [Ref. ll]. 

In a recent study Stein [Ref. 2] used the Ewing and 
Press method of solution, but solved the problem with 
the aid of modern computer algorithms. The bulk of the 
study was that of the propagation speeds of the various 
modes. His data taken from Arctic and freshwater lake 
ice experiments corraborate with the predicted transient 


arrival times and propagation losses. 


C. WAVES IN A SOLID 


There are several wave types that may propagate ina 


solid and others that occur at a solid-fluid interface. 
A study of the applicability of each type found in the 
literature was done and is included here to give the 
reader a sense of wave types that are relevant to the 
problem of modeling a sound source mounted on the ice 
generating a sound field in the underlying water by the 
excitation of the ice. 

Propagation of waves in solids can be distinguished 
from that in fluids by the solid'’s ability to restore 


shear. This leads to two wave equations: 
2 1 - 
Vo eae + (1l.la) 
p : : 


for the scalar displacement potential which describes 


longitudinal motion, and 
2 1 4 
7 es os (1.1) 
S a * 


for the vector displacement potential which describes 
shear wave motion. The potentials ¢(x,y,z,t) and 


> 
w(x,y,z,t) are related to the particle velocity vector 


~ o& 


+ ms 
u ® Ugl + Uy, J ~ Use Dy 


+ + 
u = Vo + Vxyp (isgs 
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An excellent depiction of some of the different 
types of acoustic waves, along with their speeds, is 
given by Ver and Holmer [Ref. 12, p. 273] and it is from 
this source that Figure 1.1 is taken. 

Longitudinal waves are waves that propagate along 
the direction of compression and dilatation in the 
plate. The thickening and thinning of the plate seen in 
Fig.l.la is a consequence of a non-zero Poisson's ratio. 

Shear waves are solutions to the additional wave 


equation that governs waves in solids: 
oa, a i 0 
a aa er 


where c,.% == ae in ia. Jab 


y = vector potential ‘of particle displacement 

in a solid. 
Two types of shear waves are found in the 
literature: SH and SV waves are the designations given 
them in texts on seismological wave propagation. SV, 


shear waves polarized vertically, denotes waves’ for 
which the shearing motion is perpendicular to the plane 
of the pletion It is this wave that is depicted in Fig. 
lb. SH, shear waves horizontally polarized, denotes 
waves for which the shearing motion is parallel to the 


plane of the plate. 


ry 
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Figure 1.1. Wave Types in Solids and Their Speeds 
(After Ver and Holmer) 


Young's modulus, which relates the stress _ to the 


« 
il 


strain produced, 

o = Poissons ratio,which relates lateral expansion to its 
longitudinal shortening due to a compressive stress, 

G = the shear modulus, which relates shearing stress to 
shearing strain produced, 

I = the moment of inertia, 

K = a torsional stiffness factor, which relates a twist 
to the shearing strain produced, © 

o = the solid's density, 

t = the solid'’s thickness, and 


w = the radian frequency of the motion. 
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Torsional waves, Fig. le, result from a twisting 
force applied to a solid. This wave can be generated in 
bars, where the y and z dimensions are similar. 
Torsional waves are not be easily excited in a plate; 
they will not be further considered in this study. 

Flexural waves, or bending waves, are a more complex 
motion than the previous ones considered, because of the 
complicated stresses set up when a plate is bent. The 
bending of a plate compresses the plate on the inside of 
the bend and expands the plate on the. outside. Due to 
this Smeeseion and expansion, a Poisson's ratio effect 
1s introduced. 

As with any wave, the energy of a flexural wave is 


composed of two parts, potential energy and kinetic 


energy. The potential energy is made up of two 
deflections: flexural deformation, and shear 
deformation. The kinetic energy also has two 


contributions, linear (or translatory) inertia, and 
rotatory inertia. 

The classical two-dimensional theory of flexural 
motions in thin elastic plates leads to the wave 


equation [Ref. 4, p. 31] 


2 
DV‘w + ph sy = 0 (1.3) 


3 
where D = , and here 


ILS, 


h plate thickness 


W the plate particle motion 

This equation relates only the flexural deformation to 
the linear (or translatory) inertia as the two opposing 
forces. The Mindlin-Timoshenko equation relates the 
linear and rotatory inertia as the driving forces to the 
flexural and transverse shear deformation as the restor- 
ing forces. (This was what was referred to earlier as 


"correcting" for rotatory inertia and transverse shear.) 


The Mindlin-Timoshenko equation is 


| (1.4) 
ie h3 62 3 
(v2 - 2; <a) (pv? - = “e) 2 t+ eee 
where 
at wi x 
12 2(1+o2) 
h = plate thickness. 


2 
The £. 27 term is the contribution of transverse shear. 


ph? 9 


? 
The i SE term is the contribution of rotatory 


e 


inertia. 


Note that if these terms are neglected, the classical 


equation (Eq. 1.3) results. 
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The Mindlin-Timoshenko theory is used when analyz- 
ing "thick" plates, i.e., plates whose thickness exceeds 
\,/20 , where A, is the shear wave length. [Ref. 13, p. 
2255.) 

Of key importance to this study is the effective 
phase speed for flexural motions as a function of 
frequency. Ross, assuming a harmonic disturbance of a 
homogeneous plate in vacuo, gives an equation for the 


flexural wave speed cy aS 


2 - 2 2 2 
(1- Tp(2)) (2) ° + (rtd (2) (82) - (A) = 0 (25) 
P > -B “p 
where 
lp = the plate shear parameter 
a 3 w 42 
Qp = the reference angular frequency, 
- £p 12 
" t 
w = the angular frequency of the flexural wave 
cp = the plate longitudinal wave speed, given in 


Pig. lsla. 


Equation 1.5 can be solved by applying the quad- 
e . s 2 
ratic equation to solve for the variable yy) ~— Be Se 
Pp 


a difficult function of the other parameters, but can be 
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represented accurately in terms of approximate solutions 
at low, intermediate and high frequencies. 
At low frequencies, the flexural wave speed is 
approximated by 
2 Zh 
2 w“Yt 
col * @! 136(I=o2) a9 
which is that given in Fig. 1.ld, and represents the 
flexural wave speed in thin elastic plates, given by the 


classical wave equation, Equation 1.3. This equation is 


valid when 


In the intermediate frequency range, the flexural 


wave speed is approximated by 
Soi oo | (1.7) 


In the high frequency range, the flexural wave phase 
velocity predicted by the Mindlin-Timoshenko theory 
tends to the Rayleigh wave velocity [Ref. 13, p. 248], 
which is a value slightly less than the shear wave 


velocity. The expression for the high frequency 


flexural wave speed is 


This equation is valid when 


_ 
l'p ro >>l. 


Heretofore we have discussed waves propagating ina 
solid. abe can also propagate along the boundary of a 
solid-fluid interface in much the same manner as do sur- 
face waves on the ocean. These waves are Rayleigh waves 
(for a half-solid), Lamb waves (for a plate), Stonely 
waves, Love es, and Crary waves. 

Rayleigh waves ". . .erepresent elastic perturbations 
propagating near the tres boundary of a solid and decay=- 
ing with depth" [Viktorov, Ref. 14, p. lJ]. Figure 1.2 
is a depiction of the instantaneous deformation associ- 


ated with a Rayleigh wave. 





Figure 1.2 Depiction of a Rayleigh Wave (after 
Viktorov) 
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Rayleigh waves propagate along the surface of a 
solid at a free boundary, such as an air-ice boundary or 
an air-plexiglass boundary. The wave decays exponen- 
tially in the z-direction into the solid, and if the 
decay constant is not large, the motion of the wave may 
appear at the lower (in this case, water) boundary. 
Furthermore, the Rayleigh wave is supersonic with 
respect to the fluid, therefore this motion would radi- 
ate sound into fluid (see Chapter III for an explanation 
on the importance of relative sound speeds to radiation 
into the fluid.) 

The discussion of Rayleigh waves 1s included 
because, as mentioned above, in the high frequency limit 
flexural waves degenerate to Rayleigh surface waves for 
a plate in vacuo. However, Stein shows that at low fre- 
quencies there is no coupling of the Rayleigh wave with 
the flexural wave there would be for a plate in a 
vacuo!. Moreover, the types of sources being considered 
would not generate Rayleigh waves, hence they will not 
be further considered. 

Lamb waves "“. . . refer to elastic perturbations 


propagating in a solid plate (or layer) with free 


lThere is a growing debate among the cognoscenti 
concerning the importance of the Rayleigh wave in Arctic 
ice ([Ref. 15]. 
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boundries, for which displacements occur both in the 
direction of wave propagation and perpendicular to the 
plane of the plate.” [Ref. 13, p. 67]. They are also 
called the “normal modes" of a plate, and can be sym- 
metrical (s) or antisymmetrical (a) relative to the 
center plane of the plate. Note that the zero-order 
symmetric and antisymmetric modes are the longitudinal 
and flexural waves, respectively. Figure 1.3 depicts 
the zero-order Lamb waves. These two modes are of key 


importance to this study. 


Figure 1.3 Lamb Wave Propagation (after vViktorov) 


Another wave that can exist at an interface between 
two media is a Stonely wave. According to _ Ewing, 
Jardet.zky and Press [Ref. 11, p. 112] they are modified 
Rayleigh waves that exist. when the shear wave velocities 
in the two media are nearly equal. Stein (Ref. 2, " 
69] states that they can also exist along the boundary 


between an elastic medium and a high density, 
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acoustically slow fluid or elastic medium, and that the 
Stonely wave is the high frequency asymptote of the 
first anti-symmetric mode along the ice-water interface. 
That is, they are Rayleigh waves along the solid-fluid 
iene fae. However, absorption in Artic ice is high at 
the frequencies for which Stonely waves exist [Ref. 2, 
p- 190] and therefore the Stonely wave 1s not important 
to this Study. 

Two other wave types concerning propagation of sound 
in ice are described in the literature. Love waves and 


Crary waves. Both are shear waves. Love waves are SH 


waves that are generated due to discontinuities in the.’ 


ice. Stein [Ref. 2 p. 100] notes that "“. . . unless 
there is a discontinuity ... . they will not radiate 
nor can they couple to other modes." Discontinuities in 


the ice are not considered in this study, therefore Love 
waves will not be considered. 

Crary waves are atype of SV wave discovered by 
Crary in seismic experiments on a floating ice island. 
They were generated after an ice “event", a small ex- 
plosive charge mounted in the ice. According to Ewing, 
Jardetzky and Press [Ref. ll, p. 299], this —- has the 
following characteristics: 

"l. phase velocity is near the speed of compression- 


al waves in the plate; 
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2. travel time is intermediate between that for P 
(compressional) waves and that for SH waves; 

3. principal [amplitude] is [observed] on longitud- 
tudinal horizontal seismograph; amplitudes are 
much smaller on vertical and practically absent 
on transverse horizontal; 

4. frequency is almost constant, increasing very 
slightly with time." 

Stein [Ref. 2, p. 93] again discusses these waves, 
but has found the attenuation to be at least 0.1 dB/m 
in the horizontal for the Arctic. Thus, these waves 
will not be observed at appreciable distances. Further, 
pecititire of their small vertical amplitude and velocity 
less than that of the fluid, (characteristic 3 of above) 
they do not impart acoustic energy to the underlying 
fluid. For these reasons, Crary waves will not be 
further considered. 

So far in this chapter, the types of waves that pro- 
pagate ina solid, namely Arctic ice, have been dis- 
cussed without regard to their origin. That is, they 


have been .assumed to be plane waves, far from their 


source. As these waves propagate, they may be coupled 
into the surrounding medium, creating acoustic 
radiation. 


Additional radiation may be generated when the waves 


are excited by a localized force. This is a result of 


: Zul, 


the flexural motion having a spatial gradient of dis- 
placement in the ice, which plane waves do not. Stein 
[Ref. 2 p- 96] has named this mode of acoustic radiation 
the “acoustic mode." 

Figure 1.4 is a depiction of the first antisymmetric 
mode, the flexural wave, ina solid. At this point, far 
from the source, it can be seen why the flexural wave 
does not radiate?; it travels subsonically with respect 
to the water and there is cancellation of radiation from 
adjacent positive and negative pressure peaks. (In 
actuality, because of geometric spreading, there is in- 
complete cancellation and some radiation does exist.) 
However, when these waves are excited by a localized 


~ 


force (Fig. 1.5) there is incomplete cancellation and 


-7Z 
Figure 1.4 Flexural Wave Propagation 


2A flexural wave below coincidence will not radiate. 
The coincidence frequency is that frequency above which 
the plate vibrations can couple to the fluid. See 
Chapter 3 for a further discussion. 
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radiation occurs. Note that this acoustic mode is still 
a flexural wave. A further discussion of this mode found 


in Chapter III. 
“Z 
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Figure 1.5 Flexural Wave Propagation From a Point 
Source (after Stein) 


D. SUMMARY 

In this chapter, the many types of waves found ina 
solid (ice) have been discussed. As was’7 indicated 
earlier, only a few have i importance in the Artic, 
and therefore modeling in the lab should be attempted 
only for these. Stein [Ref. 2 p. 107] has concluded 
that the three modes of propagation of importance are 
the first symmetric and antisymmetric - modes, and the 
acoustic aie All other aforementioned waves’ are 
neglected because either they are highly attenuated in 
the Arctic at the frequencies of concern or there isa 


lack of a source with which to generate them. 
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II. A LLOYD'S MIRROR EXPERIMENT 


A. INTRODUCTION 

In this chapter, an experiment will be discussed 
that tested the following hypothesis: a plate intro- 
duced at an air-water interface will influence the 
acoustic field in the water half-space. More simply, 
the reflection of a spherical wave from an Laas 
boundary develops an interference pattern with the 


original incident wave. It is hypothesized that the 


plate will effect this LI@YdJcemi tror jpattern. 


B. THE METHOD OF IMAGES 

The method of images is a well-known technique for 
predicting the radiation pattern of an acoustic source 
near a fluid interface. For continuity, a brief de- 
scription of the theory of the acoustic Lloyd's mirror 
will be presented. The following description is based 
on that given by Sanders and Coppens [Ref. 16]. 

If a source of spherical waves is placed on the z- 
axis a distance d below a pressure-release boundary, as 
shown in Fig. 2.1, the acoustic radiation from this 


source in the absence of the boundary can be specified 


as Ps = = ej(wt-kr_) (2.@R) 


where r. = y(h-a) 2 + R2 
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As a result of the boundary, there must be an addi- 


tional wave present to describe the reflection: 


pr = - a ej (wt-kr,) (2.02) 
where ry = q (nea)? + R2 
The total field is, then: 6203") 
a a 1 -jkreo 1 -jkr jut 
Ps t Pr =P =A (z-@ : =r, + Je 


The acoustic field will show strong interference effects 
because the two waves p, and p, must travel along 
different paths to reach the receiver. 


The geometry is simplified if the assumption 
r>>dcos @ (2.04) 
is made. Then, r, and r_ can be replaced by 


Goa CPAT (2.05) 
r. ™ r-Ar 
Ar » dsiné 


with 2.04, 2.03 becomes > 
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Figure 2.1 Lloyd's Mirror Geometry 


Thus when a simple harmonic acoustic wave in the the 
water is incident on a plane, perfectly reflecting pres- 
sure release boundary, the reflected wave is a mirror 
image of the incident wave which leads to an interfer- 
ence pattern in the medium below the boundary. 

This interference pattern may be observed by dif- 
ferent methods. The receiver can be placed at a _ fixed 
horizontal range from the source and moved downward to 
show the nodes and antinodes from the pattern. It might 
also be placed at a fixed slant range from the source 


and swept through an arc about the source. Or, it might 


~ 


be held at a constant depth and moved away from the 
source in horizontal range. Lastly, and want was done 
in the experiment, the source and receiver can be held 
ina fixed position and the frequency of the source 
swept, observing the response of the receiver. The 
interference of the incident and reflected waves produce 
maxima and minima in the caduenty response of the 


receiver. 


Soest EFFECT OF THE PLATE 

The effect of the plate? on the incident wave can be 
best seen by reverting to simple ray theory. Consider 
the following oe an incident wave in water, re- 
presented by ray pj, impinges ona plate of different 
characteristic impedance Fa , as shown in EP RGiene2 «2 0 ine 
ray (now p,) is refracted as it traverses into the 
plate. It reflects from the top of the plate at z =H 
(now pp). traverses the plate and then is refracted back 
mrt © the water aS Pr > Note that while the wave is in 
the plate, it is considered as a compressional (longitu- 


dinal) wave. 


3An infinite plate forced by an incident sound wave 
at any frequency, even below its coincidence frequency, 
will radiate sound, although a bending wave of the plate 
itself, below coincidence, will not. [Ref. 12, p. 280] 
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Figure 2.2 Ray Traversing Through a Plate 


The four waves may be represented as (considering only 


two dimensions): 


pi = Pi gas amy re ow (2.1) 
= j(wt-k sine,x-k cos6,Z) ; 
Pa = Pa @ 2 a 2 a (2.2) 
= }(wt-k sin6,xt+k cos 6,Z) , 
Pb P e J . b > b (2.3) 
is J % 
Pr = Py ej(wt e sin watee,2) (2.4). 
The boundary conditions are: 
at z = 0, continuity of pressure requires 
pi (x,o) + pr(x,o) = Palx,0) ? Pp(*,0) (2.5) 
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Also at 2z=0, the condition of continuity of normal 
velocity, requires (taking into account direction of 
propagation) 

uj (x,0)cos6, - u,(x,0) cos 6, 

(2.6) 
= u,(x,0)cos6,-u,(x,0) Cos 6} 

The boundary at z =H is pressure release, therefore 

Pa(x-H) + Pplx,H) = 0 (2.7) 


For each plane wave, we have 


2 = +r (28) 
Xu 
With Snell's Law, sin Sra sa ; (2.9) 
1 2 


and since the angle of incidence equals the angle of 


reflection, we have 
05 = oy = 0, and O5 = Opn = 05 (2.10) 
If each value of u ‘in 2.6 is replaced with the 


appropriate value of p/r , and from 2.10, subscription 


are changed, then: 


SEAT [pi (+0) -pe(X,0) ] = — [pa(x,0)-pp(x,o) ]'2°4?) 


Or 


x 9 
pilx,o)-py(x,0) = oar (p,(x,0)-pp(x,0)) (2.13) 


Apply 2.1 through 2.4 to 2.5 , and drop the eJ4t , 


-jkK sin6é x - -jk siné x 
(PitPp)e 2 1 = (Pathe 22 (2.14) 


and do the same to 2.13, 


(Py-P_)e JK181NGb1% S ae 2 (P_- Pe” Jk 28ine28 (9, eS) 
Z 


Apply the boundary condition at z=H (Equation 2.7), 


and 
-jk sine x 
dividing through by e 2 2 gives 


\ 


Algebraic manipulation of 2.14 - 2.16 then leads to the 


reflection coefficient R 


P (ie 21 8%) - 01 ta Se e7 J2k2Cc0s86 2H 

2 ome o 5 Ee a 

= Pj (1 + z1 SoS". ) a4 . <i, See) =j2k, cos 65H 
cosé a cosé 


i (2.17) 


where ?; and e are related by Snell's Law (.9), and 


cOSsQ@5 = l1- (=2)2 sine, 
l 


(2. a) 
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D. VALIDATION OF EQUATION 2.17 

In the limit as H+0O , Equation 2.17 reduces’ to 
R =-l1 . This is to be expected, because when H+0O the r 
layer disappears and the incident wave is * reflected 
off a purely pressure-release boundary, where the re- 
flected wave is 180° out of phase with the incident, and 
iWence R= -l . 

In the limit as — , that is, normal incidence, 


Equation 2.17 reduces to: 


—- 70) Cosk2H + jrogsink2H 
RL ri cosk 5H + qe Sink 5H (2 19) 
Comparing Equation 2.19 with Equation 6.12 of. 


Kinsler et al [Ref. 17] 


r ~ft r : 
es ee et (2.20) 


x (1+ =1)cosh 2H + i + [1)sinh 2H 


where a is characteristic impedance of the fluid above 
the place. In the limit — », that is, assuming a 
vacuum: above the plate, 2.20 reduces to 2.19. 

Finally, comparison of 2.12 with Brekhovskikh [Ref. 


18] yields additional validation. His coordinate system 





which reduces to [Ref. 18, Equation 5.18] 


ZJaod 
ie eau) 9 CE (2.24) 
eee ey 2S 
Bas 12e 
_ 24-22 _ 29-23 
where V}92 7 a7! V53 Ze (2.253) 
1 2 


Assuming a = 0 (in vacuo), then Ps = -l, and .27 


reduces to 


Z Z Zjaod 
=a —- ze 
a (2.26) 
™ Z Z ZjJaod 
(1 + 34) - (1 - F4)e7/% 
ae 2 
changing to previous notation and noting that d = <-H , 


equation 2.26 becomes 2.17: 


aoe 251 a) rt + Zi ogee a7 ZIK2C0S6 2H 


CT, cosé ry, cose 
(1 pa gents, =e Ei S282 a7 ZIK2¢00Ss 6 oH 
ae ee a oe 


E.- VALIDATION OF THE PLANE WAVE ASSUMPTION 

Thus far, the reflection coefficient of an oblique 
incident plane wave off a finite layer has been calcu- 
lated and validated against previously known theory and 


limiting cases. However, one question still needs to be 


addressed before this equation is put to use. That is, 
the validity of the plane wave assumption made above, 
Equation 2.8, needs to be established. 

The solution for a point source near a boundary is 
addressed in two works. Officer's [Ref. 19] solution is 
in the form of three integrals for each of the three 
regions. However, the integrals are difficult to 
compute. Officer does give ". . . the condition under 
which the plane wave reflection coefficient may be taken 


as a valid approximation for the reflection of an acci- 


dent spherical wave. This condition is that 


2 S i oe (2.27) 
w Ry (Sy - sin20]3/2 
2 


where R) is the distance from point of reflection to the 
receiver. 

The second work that addresses this problem is 
Brekhovskikh. From Brekhovskikh, the geometric optics 
approximation (in which a spherical wave is relected off 
a layer with the same reflection coefficiient as a plane 
wave) is valid when the distances of the source and the 
receiver from the boundary are large compared with the 
wavelength. [Ref. 18 p. 256] In the experiment run, 


this assumption is valid above f = 50kHz. 
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Lastly, for completeness, the extent of the near 
field should be calculated. From Ziomek [Ref. 20], the 
maximum radial extent r of the near field of the source 
is given by: 

x > aR (2.28) 

r 
where 
R is the maximum radial dimension of the transmit 
transducer. 

For the experiment done here, the following values 


were used in equations 2.27 and 2.28: 


cy = 1481 m/s 
co = 1992 m/s 
6) = 42.27° 
R; = 0.5 m 


Equation 2.27 requires that f >> 14.8 kHz for the plane 
wave assumption to be valid. From 2.28, for the radia- 
tion from the transmitting LC-10 to be in the far-field 
prior to reaching the plate requires f < 1.01 MHz. 
Therefore for the frequency range of interest (50-100 


kHz) both conditions are satisfied. 


FF. APPLICATION OF EQUATION 2.17 
Upon close inspection of Equation 2.17, it can be 
seen that the reflection coefficient is a phasor of unit 


magnitude with a phase delay that is a function of 


frequency. This is intuitively true because of the con- 
ditions imposed on the problem: the boundary condition 
at the top of the plate, Z =H , is a pressure release, 
perfectly reflecting plane; the medium is assumed to be 
lossless. Since it is perfectly reflecting, there is no 
refracted wave at the upper boundary and there is no 
loss of signal. If the medium is lossless there is no 
loss of signal as the wave travels through the plate. 
Thus, the reflection is of unit magnitude with a phase 
angle, that is R = l‘eJT, 


The phase angle +t of R Is; from Equaerten 2317, 


| 1+ Sl SOS%2)sin[2k jHcos 6, | 
tT = tan 
_ ZX, COS85)_ XZ) cos® 
(1 bs scene r Sosg )COSl 2k 2Hcos 6 2] 
(2.29) 
_ ©, cos9 . 
ae (1 - Sb SSa57)sin(2k2Hcos 62) 
(1 + ZL-S08*2 )cos (2k pHcos 82) 


In the experiment to be described, the following para- 


meters were used in Equation 2.29: 


= = kK = 6 Pas 
r) = p,c, = (998 —$] [1481 m/s] = 1.48 x 10 — 
= = k = 6 Pa’s 
ro = poc, = [1200 ~4] [1992 m/s] = 2.39 x 106 =“ 


H = 3.18 x 1073 m 


From the geometry of the setup, 


drevr - 5m 
l drevr + Isource -55m 


where droyr and Qcource 22e the receiver and source 
depths, respectively. A plot of the phase angle caused 
by the plate, Equation 2.29, as a function of frequency 
with the parameters shown in Figure 2.4 above is in 
Figure 2.4. Note the conventional of a negative angle 
producing a phase delay in the signal. 

In the frequency range of interest, f = 0 to 100 
KHz, +t decreases from O° at f£ = OHZ, to -104° at 


£ = 100,000 kHz. 
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Figure 2.4 Phase Delay Caused By the Plate 
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The phase delay caused by the plate will be mani- 
fested by a shift in the frequency of the nodes and 
antinodes of the interference pattern. It can be shown 


that the frequency shift Af caused by the phase delay rt 


1s; 
Cot 
where 
Ar = difference in path lengths between the direct 
and reflected signals 
t = phase delay caused by the plate, as obtained 


from Equation 2.29, 


Co speed of sound in the fluid. 


G. DESCRIPTION OF EXPERIMENT 

Two small LC-10 transducers were placed in the 
water, one acting as a source, the other aS a receiver, 
and were connected to a Hewlett Packard 4162A frequency 
analyzer as shown in Figure 2.5. A periodic chirp sig- 
nal of 0-100 kHz was applied to the source transducer, 
and the received response was measured and plotted. The 
LC-10 transducers were assumed to have a flat response 
throughout. the frequency range. While the response 
curves of the LC-10 show the response to be fairly flat 
up to 100-kHz, the assumption is really not tested, 


Since the response of the transducers themselves 


is not being tested, rather the difference in response 
with and without the pressure of the plate. 

. The response of the receiving transducer showed the 
maxima and minima in the pressure at its location due to 
the boundary. The response was measured for the water- 
air boundary, and without changing the location of the 
two transducers, the plate was introduced and the 


response was remeasured. 


H. EXPERIMENTAL RESULTS 

The results of the experiment can be seen from 
Figures 2.6 and 2.7. Figure 2.6 is a plot of the fre- 
quency response of the receiving hydrophone with fixed 
geometry with no plate present; the classic Lloyd's 
mirror experiment. Note the frequency pattern of the 
aes and antinodes. When the plate was placed in posi- 
tion without moving source or receiver, the positions of 
these nodes and antinodes shift to the right, higher in 
frequency, as can be seen in Figure 2.7. This was pre- 
dicted as the manifestation of the phase angle caused by 
the plate. The change in the frequency of the nodes 
increases linearly with frequency, corresponding to the 


increased phase angle. 
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Figure 2.5. Lloyd's Mirror Experiment set-up 


Table 2.1 shows a comparison of the predicted nodal 
shift to the nodal shift obtained experimentally. ™ was 
obtained from Figure 2.4. The predicted frequency shift 


was obtained from Equation 2.30. 





TABLE 2.1 
tnode T Af predicted Afaxp. 
37 kHz -42° eciiz, iao2 KHZ + .13 
56 -62 digs ZO + .43 
5 -81 1.9 2.24 fe a 
Jo -100 2 2-50 tine OO 


The data of Table 2.1 are shown in Figure 2.8. The 
size of the error bars of the experimental data are due to 
the inaccuracies in reading the frequencies of the nodes. 
Note, however, that in all cases the experimental data 
agreed with prerdicted values of the frequency shift with- 
in the accuracies of measurement, so on a qualitative 
level, at least, the theoretical predictions of the nodal 
shifts are supported. 

In further comparison of Figures 2.6 and 2.7, it can 
be seen that although the dictern has been shifted in 
frequency, there 1s no change in the amplitudes of the 
maxima and only a small change in some of the minima of 
the frequency response pattern when the plate is 
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Figure 2,6, 


free boundary 


Frequency response, 


introduced . This suggests that there is insignificant 
change in the amplitude of the reflected signal, as pre- 
dicted by the refllect¥on coefficient KR Maving a4 ees. 
of one. From this, it is inferred that there is insignif- 
icant translation of the incident wave into plate wave 
modes other than longitudinal. The ray-~-method used in de- 
riving the reflection coefficient R assumes that the com- 
pressional wave incident at the plate boundary remains a 
comressional (or longitudinal) wave in the plate. 

A change in the mode of the wave from compressional 
to, for instance, flexural or shear, would involve a loss 
of energy in the reflected wave, which would be seen asa 
Change in the amplitude of the maxima and minima of the 
frequency response pattern. Since very little change was 
observed, it appears that the incident wave remained a 
compressional wave as it traversed the plate and re- 
entered the fluid below. 

Thus, it has been demonstrated that the effect of 
placing a plate into a fixed-geometry a tioa~o fined 
problem is to introduce only a phase delay into the re- 
flected signal, causing a related upward shift in the fre- 
quencies of maxima and minima in the pressure for the 


fixed source and receiver positions. . 
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III. RADIATION FROM A POINT DRIVEN 
FLUID~-LOADED PLATE 

A. INTRODUCTION 

Chapter II treated the effects on the sound field 
from a source in the water when a plexiglass plate is 
Placed on the surface. This chapter will deal with the 
acoustic radiation from a plate driven by a_= source 
attached to it. When the source is driving the plate, it 
excites the "acoustic mode" of sound propagation in the 
plate, as described in Chapter I. From the previous 
discussion, the acoustic mode leads to radiation away 
from the plate caused by the spatial gradient of the 
amplitude of the flexural wave excited by the source. 

The principal references for the theoretical de- 
scription of the problem are Feit [Ref. 7] and Ross 
[Ref. 21]. A computer program was written to compute 
the radiation pattern from the driven plate at various 
frequencies. An experiment was conducted to obtain 
actual radiation patterns from a plexiglass plate placed 


on the water. 


B. THEORETICAL DEVELOPMENT 
As was seen in the background section of Chapter 
One, the theory concerning the point-driven fluid-loaded 


plate is well known, and its full development will not 
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be given here. However, for continuity, key points 
leading to an understanding of the solution will be 
summarized. 
1. Radiation from a Propagating Flexural Wave 
A flexural wave propagates in a plate with wave 


number Kr where 


When this plate is in contact with a fluid, vibrations 
of the plate will be imparted to the fluid, producing 


waves in the fluid of wave number Ko where 


+> +> 
ko = k, + k, (3.2) 


> - 
where k, and k, are the wave numbers perpendicular and 


‘parallel to the plane of the plate, respectively. 





be 


ot —) Ge = te +a ee 





Figure 3.1 Sound Radiation from a Flexural 
Wave of a Plate 
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The length of this vector is given by 
ko* = koa (43539) 
the component parallel to the plate is 


and from Snell's Law this must be equal to Ke - Trae 


leads to [Ref. 20, p. 163] 


From the boundary conditions at the interface, the 


normal components of particle velocity must be equal, 


Vv = u (3.6) 
te Ly 


where ve refers to the perpendicular component of the 
plate particle velocity and ~.. is the fluid particle 
velocity perpendicular to the plate. 

From these conditions, it can be shown that [Ref. 21] 
p- 163] the acoustic pressure Pp on the plate due to 


the fluid is: 


-JWwP—olo ———— ¥ (35 7) 
yo - Kg 


where w is the normal displacement of the flexural 


e 


wave. Equation 3.7 gives the acoustic pressure p on 
the plate attributable to fluid loading. Note that p is 


dependent on the relative values of kp and kg . 
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This leads to the important effect of coincidence. 
When k¢>Kk, , the expression for the fluid pressure is 
real, which represents mass~loading [Ref. 21 p. 163] and 
acoustic power is not radiated into the fluid. When 
ko = Kg . the condition of coincidence is reached. The 

expression for the radiated acoustic pressure from the 
flexural wave would seemingly predict infinite pressure. 
(This cannot, however, occur, because fluid loading pre- 
vents the two speeds from being exactly equal [Ref. 21 
p- 163]) When kg<k, , the expression for the fluid pres- 
Sure in Equation 3.7 is imaginary, so the effect is that 
of a radiation resistance, [Ref. 21 p. 163] and the pro- 
pagating wave Pie. power away from the plate. 

The coincidence frequency is obtained by equating 
the wave numbers of the flexural wave in the plate and 
the compressional wave in the fluid and solving for fre- 
quency wo - But because, the flexural wave speed varies 
with frequency and formulae for the flexural wave speed 
depend on whether. thin plate or thick plate approxima- 
tions are used, there are different forms for the solu- 
tion for the coincidence frequency. 

| As there were three approximations for the flexural 
wave speed (Equations 1.6 ~ 1.8), there are three 


formulae for the coincidence frequency. In the low 


a2 


— 


frequency limit, or thin-plate theory, the coincidence 
frequency is given by: 
2 
wey = So (3.8a) 
Pp 
In the intermediate frequency range, the coincidence 
frequency cannot be easily solved for implicitly, but is 


given by the solution of: 











Woz eye Wot 
wei = (1 + (G—)2)  - 2k (3.8b) 
[SS 2 
Q Wr GC 
p C1,2 Pp 


In the high frequency range, the coincidence frequency 


Can be given by: 


2 
wen * Sp 120 __ SF (3.8c) 
- t ea" 
Cp 7 


Note here that if Cp <¥2 Sy > ‘Equation 3.8c is 


negative, which implies that coincidence cannot be 
reached. Recall that the high frequency asymptote of 
the flexural wave speed is the Rayleigh wave speed, a 
value somewhat less than the shear wave speed. If this — 
speed is less than the speed of sound in the fluid, the 


condition of ky = Ke, or coincidence, cannot be reached. 
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2. Radiation from a Point Source 

Heretofore, we have discussed the flexural wave 
itself radiating power from the plate. The conclusion 
that no power is radiated below coincidence applies only 
to infinite plane bending waves on an infinite plate. 
If excited by a localized force even below coincidence, 
sound will be radiated because the waves in the plate 
have divergence [Ref. 20 p. 169]. 

The pressure at a field point (R, 6, 9) froma 
point-driven fluid-loaded plate can be, in a derivation 


by Feit [Ref. 7, p. 1492], eventually expressed in the 


form - (3.9) 
P(R, 6, 8)= Ea Ka st) (x oRsingsinw) eJ*oRS1N$COSWy, (1) 5 i nwaw 
C 
where 
p(R,¢,9) = acoustic pressure in spherical 
coordinates with the source at the 
origin, 


R = the slant distance from the source, 
¢ = the azimuthal angle, 


§ = the circumferential angle measured in 
the plane of the plate, 


w = the variable of integration, 
Ky = the wave number in the fluid, 


C = contour of integration shown in 
Figure 3.2, 


Fo = the amplitude of the source. 


ie 
- = the zero-order Hankel functron of the 
Oo first kind, and 


V(w) = a complicated formula that relates the 
vibrating field of the plate to the water 
through the Mindlin-Timoshenko thick 
plate equation. V(w) contains the poles 
of the integral in Equation 3.10. 

The solution of the integral was obtained by 

Brekhovskikh (Ref. 18, pp. 242-258] and is given by Feit 


[Ref. 7, Equation 24]. Ross, however, gives a more use- 


ful form of Feit's solution: 


= TkoFo COs 
p(R, o)= SoMa o - (3216p 
1-3 cos a sin2¢-(c sin*o-._(c 
W Cc ee a ee 9 
8 1+() (eo) [p(sin29 (en) 


where Ip is the plate shear parameter, given by 


~ 2.65(1+ 3 o + = a?) (3,28) 


'p 
and 8 is the fluid loading factor, given by Equation 
3.13, which is discussed fully in the next section. 

The bracketed term in Equation 3.10 along with 
the cos¢ can be described as a directivity function 
H(w,¢)- This directivity function was encoded in 


FORTRAN and the program (given in Appendix A) was’ used 


to generate plots of theoretical directivity patterns 
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at several different frequencies. From these figures 


(Figures 3.5-3.11) the coincidence effect can be readily 
irvea. The figures were based on thin plate theory, 
and cover frequencies from below to above coincidence. 
They represent a hypothetical plate and are presented 
for a discussion on the effects of coincidence and fluid 
loading. They will be discussed in section El of this 


chapter. 





Figure 3.2 Integration Path of Integral of Equation 3.9 
(from Feit) | 
C. THE EFFECT OF FLUID LOADING 
Although Equation 3.10 is valid for all ranges of 
plate and fluid parameters and frequencies, the expres- 
sion may be simplified greatly and further insight may 


be obtained by considering the effect of fluid loading. 


=) 


The fluid loading factor, 8, 1s the ratio of the 
specific acoustic impedence of the fluid to the mass 
reactance per unit are of the plate. Ross [Ref. 20 p. 


166] defines it as 


O S 

B = ——- (ani) 
WPpD 

where 

Po = density of the fluid 

Co = speed of sound in the fluid 

w = radian frequency of the flexural wave 

Pp = density of the plate 

t = thickness of the plate 


Applying the formula for the thin plate coincidence 


frequency, Equation 3.8a, the fluid loading factor 


becomes 
V12 Pp So 8 
where 
Cp = the plate longitudinal wave speed, given in 
Figure l.la, 
Wo = the coincidence frequency. 


A fluid loading factor much less than 1 corresponds 
te light fluid loading, which occurs, for example, ina 
metal plate in air. A fluid loading. factor much greater 
than 1 corresponds to heavy fluid loading found, for 


example, a plexiglass plate on water. The effect of 
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fluid loading can best be seen by considering four 
general cases: frequency below coincidence, light and 
heavy fluid loading; frequency above coincidence, light 
and heavy loading. 

When w<wo and B<<l, a case of light fluid loading, 
the radiation pattern is that of a small baffled piston, 
or monopole-type radiation, non-directive. For this 


case, Equation 3.10 reduces to [Ref. 7 p. 1492] 


90F 
ilea, d, 0) ar eJkoR (3.14) 


Ross [Ref. 20 p. 170] states that the radius of this 
piston is the sum of the radius of the surface over 
Witch the force is applied and a quarter of the flexural 
wavelength. Ver and Holmer [Ref. 12 p. 278] give a 
slightly different form. They state that for a plate 
driven by point source, the equivalent radius of the 
baffled piston is 0.286: Aq » where Aj, is the flexural 
wavelength at coincidence. 

When, w<wo and g>>1 , heavy fluid loading, the point 
force acting on a thin plate radiates sound as through 
the plate were not there, representing dipole-type 
radiation, and the radiation pattern is cosinusoidal. 
Equation 3.10 reduces to [Ref. 7 p. 1492]: 


~IKoFo @JkKOR 


aR COS 6 (3.15) 


p(R, b ) ~ 


When w>w, and g<<l , the radiation pattern is still co- 
Sinusoidal [Ref. 21 Gon sigag). 

When w>wo and 8?>1 , Equation 3.10 can be reduced 
only slightly, and gives the interesting radiation pat- 
terns corresponding to the flexural wave radiation and 


radiation from the driving piston. 


D. THE EFFECT OF STRUCTURAL ABSORPTION 

Above coincidence, then, the absorption of the flex- 
ural wave as it propagates through the plate will effect 
the radiation pattern. Feit [Ref. 7, p. 1494] states 
this and states that (without proof) the far field pres- 
sure maximum 1s. reduced by 5 dB for a loss’ factor, 
n=0.1. (See Equation 3.18 for the definition of the loss 
factor n'.) Ungar [Ref. 22, p. 453] lists a loss factor 
typical of plexiglass of n = 0.02 to 0.04 in the audio 
frequency range. Therefore, it 1s expected that the de- 
crease in o_o peak pressure amplitude due to the 
absorption of the flexural wave will be less than 5 adB. 

An experiment was conducted to quantify the loss 
factor for the plate used in the experiment. This was 
done also to investigate the effects of the refléction 


of the flexural wave from the edge of the plate, l.e., 


*Although the effect is the same, the term absorption 
1s perhaps more correct than damping as the term damping 
implies additional materials applied to the plate. 
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at “88 ious tes, POlnts 
“OuPling between 


crucial] and 
end 


patterns were induced. (Other samples, for instance 
of metal, did produce Schladni patterns using this 
technique. ) 

Qualitatively, the effect of the plexiglass loss 
factor n was very apparent. Refer to Figure 3.3 which 
is a plot of flexural wave attenuation, after correction 
for cylindrical spreading. As the accelerometer was 
moved farther away from the source, the corrected ampli- 
tude of the flexural wave steadily decreased with range, 
dropping 40 dB by the end of the plate for the frequen- 
cles of interest. No reflections of the pulse off the 
edges were observed. 

Ver and Homer [Ref. 12, p. 276] give the losses in 


power radiated from the plate due to absorption to be 


ec TORE R/2 
where kf = flexural wave number (3.18) 
R = range from rource 


r loss factor of the plate material 


using this, with the 100 kHz absorption curve (the only 
one straight enough to be trusted) yields n ~ 0.03 
which is consistent with the values of Ungar [Ref. 21, 
p.. 453]. 

According to Junger and Feit [Ref. 13, p. 249], the 
effects of structural damping on the far-field pressure 


pattern are accounted for, for small loss factors, by 
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multiplying ¢c, in Equation -3.10 by (l-jn)}/2 
~ (1-1/2 jn), or we by (1+ 1/2 jn). They further state 
that nmeniied plays no significant effect. on the maximum 
pressure of the radiation field below coincidence. This 
is also intuitively true because below coincidence the 
propagating flexural wave does not radiate and therefore 
its absorption has no effect. Above coincidence, they 
state (as Feit did in his earlier work) that the effect 
is to reduce the far-field pressure maximum. Figure 
3.4, taken from Feit [Ref. 7, p. 1494], show the effect 
of structural damping on the far-field of a point- 


excited thin plate. 
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Figure 3.4, The Effect of Damping on the 
Far-Field Pressure (from Feit) 
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. E. DISCUSSION OF THEORETICAL PLOTS 
1. The General Effects of Increasing Frequency 

Figures 3.5 to 3.11 were generated from the pro- 
gram given in Appendix A. They show the radiation 
directivity pattern H($), normalized to the on-axis 
pressure, for a frequency range from below coincidence 
to above coincidence. 

The plate parameters used in these figures were 
Bmose of a hypothetical plate and fluld combination. 
They do, however, show quite well the effects of in- 
creasing frequency on the radiation pattern in the 
fluid. Thin-plate theory was assumed, and the ratio of 
w/wo Was the parameter that was varied. 

Figures 3.5 to 3.7 show the directivity patterns 
below coincidence. Far below coincidence,Figure 3.5,the 
radiation from the plate is cosinusoidal, with maximum 
radiation at ¢ = 0° and minimum radiation at 6 = 90°. 
It was shown in Section Bl of this chapter that be low 
coincidence an infinite plane propagating flexural wave 
will not radiate sound. Tt was shown in Section C that 
the radiation pattern should be that of a low frequency 
small radius piston mounted on an infinite baffle. This 
ils revealed in these first three figures. Note that as 
the coincidence frequency is approached, the pattern 


changes into a wider (Figures 3.7 and 3.8) pattern. 


As the frequency of the source becomes larger 
than the coincidence frequency, the phase speed of the 
propagating flexural wave exceeds that of the fluid, and 
therefore the propagating flexural wave can radiate 
energy away from the plate. This is seen in Figures 3.9 
and 3.10. 

| Well above coincidence, Figure 3.11, a third 
lobe is predicted. This, according to Feit, is_ the 
". . »« nonpropagating near-field deflection of the 
plate .... itself becoming a propagating wave." 
[Ref. 7, p. 1493] 

Feit, in a later work with Junger [Ref. 13, pb. 
248], state that thick plate theory yields larger values 
for the angles at which the propagating flexural wave 
radiates and smaller values of peak pressure. That is, 
the lobes in the pattern will be at a greater angle and 
have a smaller value. Figure 3.11, from Feit and 
Junger, show a comparison of the size of the one 
lobe and the angle at which it appears for thin-plate 


theory and thick-plate theory. The figure neglects the 


effect of structural damping. 
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2. The Effect of the Fluid Loading Factor, 8 


The plots of the previous section showed the 
general trend in the radiation pattern when frequency is 
increased for a given set of thin-plate parameters. 
However, there are subtle changes that take place in the 
patterns when different parameters are used. The den- 
sity of the plate can be measured with reasonable 
accuracy, and Poisson's ratio should be somewhat similar 
among plexiglass samples. However, Young's modulus and 
the exact effect of Poisson's ratio in the wave speed 


formula may vary. Additionally, Equation 3.10 lends 
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itself better to changing certain parameters than 
others. The parameter chosen for variation was’ the 
fluid loading factor, g , for fixed ratios of w/w. 
Figures 3.13 to 3.21 are the results of these runs. The 
effect of varying the fluid loading factor g as follows. 
Increasing the fluid loading smooths the radiation 
pattern, reducing the peaks and making the coincidence 
effect less dramatic. 

Below coincidence, fluid loading should fot have 
a great effect on the radiation pattern. This is seen 


Pieregures 3.13 thru 3.15. A change in the fluid load- 


ing factor g for a frequency ratio of wo/w = 2 below co- 
incidence has a minimal effect on the radiation 
pattern. 


At setlactiaenas, the effects of fluid loading on 
the radiation pattern can be seen: for w@iatght fluid 
loading, Figure 3.16 shows a well-developed radiation 
pattern, but as g is increased (Figures 3.17 and 3.18) 
the radiation pattern is less pronounced and the off- 
axis pressure gain not as great. 

Similarly, above coincidence, (wo/w) = Q.86, as 
seen in Figures 3.19 -through 3,20 the increase in fluid 
loading causes a less dramatic radiation SaMXDn . 

The importance of g is noticeable. Even though 


in the computer runs £g was never "much less’ than one" 
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which was described as the condition for light fluid 
loading, small changes in £8 caused large changes in the 


radiation patterns. 


F. DESCRIPTION OF EXPERIMENT 

An experiment was conducted to measure the radiation 
pattern from a point-driven sheet of plexiglass. See 
Figure 3.22. A 1/8th-inch thick sheet of plexiglass of 
dimensions 87cm x 152cm was floated on the surface of a 
freshwater anechoic’ tank. A small silicon bead was 
applied to the edges of the plate to keep it afloat. In 
the center of the plate a Tonpilz longitudinal vibrator 
of the author's design5 and manufacture was bolted to 
the plate. The vibrator was driven by an HP-3314A 
function generator operating in the pulse mode . A maxl1- 
mum driving output voltage of ‘30V was applied. 

A schematic cross-section of the Tonpilz vibrator is 
given in Appendix BB. Its tail mass was made out of 
stainless steel and its head mass out of aluminum. This 
was to Mbart as large a signal as possible to the a 
iglass sheet while keeping the vibrator's natural 
resonance frequency fairly high. An admittance circle 


of susceptance vs. conductance and plots of admittance 


“For an excellent description of the design of 
Tonpilz vibrators, see Wilson [Ref. 23, Ch. 6]. 
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and susceptance vs. frequency for the vibrator are also 
given in Appendix B. They show the motional resonance 
frequency of the vibrator to be 92.4 kHz, anda Q of 
about 112. 

It is stated that although the motional resonance 
frequency of the vibrator was within the frequency range 
of interest for this | experiment, this was. not of 
key importance because a swept-—frequency mode was not 
used. However, what is of interest in the pulse mode is 
the Q of the vibrator. The Q of the vibrator showed its 
effect in the transient response of the vibrator in the 
pulse mode. The transient response decay time grew when 
the driving frequency approached the motional resonance 
frequency of the vibrator. However, in every data set 
measurement, the pulse length was AG aiaill to give maxi- 
mum Steady state response time for the given pulse 
interval. 

The head mass of the vibrator was about 1l cm in 
radius and was attached to the plate with a thin film of 
Silicon grease to ensure positive contact between the 
head mass and plate. The clamp rod of the vibrator was 
used to compress the vibrator against the plate. 

The radiation pattern from the driven plate was mea- 
sured with a Celesco LC-10 hydrophone whose output was 
sent through a Tektronics AM-502 low noise preamplifier 


to a Kikusui oscilloscope. The LC-10 was mounted on an 
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HP 3314A Kikusui COS 5060 
Function Gener ator Oscilloscope 









Tektronics 
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Accelerormeter 


Plexiglass Plate Nm 
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; LC-10 
Hydrophone. 


Ahechoic Tank 


Figure 3.22, Driven Plate Experimental set-up 
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angular jig that allowed it to be swept from 4 = -90° to 
90° at constant range from the vibrator. An Endevco 
K920 accelerometer was used at times to probe the plate 
to obtain attenuation and edge-reflection data. 

The pulse mode was used because of a problem of 
electrical pickup between the leads to the vibrator and 
the LC-10. The electrical pickup could not be removed 
even with the use of grounded coaxial cable and shielded 
connectors. 

The length of the pulse varied with frequency. The 
limiting factors for the pulse length were the pulse in- 
terval and the transient response time of the vibrator. 
The pulse interval was limited because the LC=10 was 
located at a fixed distance from the source and recep- 
tion of a pulse by the hydrophone had to be ended before 
the next began to prevent overlapping signals. The 
transient response time of the vibrator also limited the 
length of pulse used. The received pulse showed the 
transient and stead-state response of the vibrator. 
However, amplitude data could only be taken on the 
Steady-state portion of the received signal. This 
further limited the time available in the pulse with 
which to drive the plate in a steady-state condition. 

Both the pulse interval and pulse length were ad- 


justed for each data set measurement to provide the 
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maximum time available with which to drive the plate in 
a steady-state condition. 

The properties of the plexiglass used were deter- 
mined experimentally by Johnson and Denny [Ref. 23|] 
Table 3.1 lists the wave speeds for the three wave types 


they measured. Note that both the shear wave speed and 


Experimental Wave Speeds 


for Plexiglass 





Cp = 2240 m/s 

Cs = 1373 m/s | 

Ch = 700-1150 m/s | 
. | 


therefore the -high-frequency asymptote of the flexural 
wave speed are less than the speed of sound of the 


water. Therefore, coincidence cannot be reached. 


G. EXPERIMENTAL -RESULTS 
1. Discussion of Plots 

Several plots were obtained at frequencies from 
50 kHz to 245 kHz with varying degrees of success. As 
was stated earlier, the high frequency asymptote. of the 
flexural wave speed was less than the speed of sound in 
the fluid, so the coincidence effect could not be 
obtained. Still there was a change in the radiation 


pattern from the plate as frequency was increased. 
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Figure 3.23 shows the experimental radiation 
pattern when the plase was driven at f=50 kHz. There 
are several aberrations in the pattern when the deflec- 
tion angle increased above 40° . Two plots were made, 
one with the hydrophone at a range of 0.5 m and the 
other at a range of lm. This was done to see if the 
distance from the source to the vibrator had an effect 
on the pattern. It can be seen that it did not, as the 
two curves appear to have quite similar behavior. 
Still, the general cosinusoidal pattern of a point- 
driven plate will below coincidence is seen. 

In Figure 3.24, the aberrations of the previous 
pattern are more pronounced. The next sets of data, 
Figures 3.25 and 3.26 at 92.3 and 93.5 kHz, respectively 
show the development of a side lobe at an angle of about 
45° . Figure 3.26, where data were taken at angles of 
-~90° to 90°, show an asymmetry in the pattern. This was 
seen in a few plots and is thought to have been caused 
by asymmetrical excitation of the plate. This is dis- 
cussed further in section G.3. 

The next three plots, Figures 3.27 to 3.29, at 
frequencies of 100 kHz, 120 kHz, and 150 kHz _ show the 
development of lobes at about 40°. ° There are multiple 
aberrations in the pattern which may or may not have 


been additional minor lobes. Note that only the pattern 
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at £ = 150 kHz show a clean major and minor lobe 
pattern. 

In the following three sets of plots, figures 
3.30 to 3.32 at f = 170 kHz, 194 KHZ, .and 1% Kez eae 
well-developed pattern of Figure 3.29 is lost and the 
pattern became more complex, with multiple aberrations 
which may have represented small minor lobes in the pat- 
tern or interference resulting from incomplete and asym- 
metrical excitation of the plate. 

The final three data sets, figures 3.33 to 3.35 
at frequencies of 225 kHz, 227 kHz, and 246 kHz, respec- 
tively, show patterns that are very complex, and as evi- 
dence in figures 3.33 and 3.35, very asymmetrical. The 
side lobes of the pattern are just as large as _ the on- 


axis lobe. 


2. Error Analysis 


The method of data collection was as follows. 
For a given frequency, with the hydrophone on-axis (0°) 
the output voltage of the preamplifier was adjusted to 
7.0V. All readings thereafter were immediately corrected 
to dB re 7.0V. The signal-to-noise ratio (SNR) varied 
between data runs because of the output response of the 
vibrator. On the average, however, the SNR was esti- 


mated to be 30 dB at zero deflection angle. 
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Since the data were normalized to the on-axis 
reading, the error increased with decreasing sound pres- 
sure (and hence voltage) levels. Table 3.2 shows the 


resulting dB error for various levels. 


EXPERIMENTAL ERROR 


TABLE 3.2 
SPL ERROR 
OdB re 7V reel OB 
a) + 0.56 dB 
-10 + 0.84 dB 
-15 fez dB 


Angular measurements were quite precise with the 
jig that was used. Angles were precise to +0.1 degree. 
As precise as the angular measurements were, they were 
not as accurate. While every attempt was made to accu- 
rately position the hydrophone directly underneath the 
vibrator at the start of each data collection period, 
parallax and the water refractive effect hampered accu- 
rate positioning. It is estimated that the angular 
accuracy was +5 degrees. Techniques to improve the 


accuracy were abandoned because it was reasoned that the 
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the result of any inaccuracies would be an angular 
skewing of a pattern, not a change in the pattern 
itself. 

It would seem therefore, that inaccuracies or 
errors in data collection were not the cause of such 
large aberrations are asymmetries in the experimental 


plots. 


3. Discussion of Results 

Since the high frequency asymptote of the flex- 
ural wave speed was less than that of the speed of sound 
in the fluid, coincidence was not possible. Therefore, 
the large off-axis pressure maxima associated with the 
driven plates above coincidence were not obtained. This 
may, however, be proper modelling of the ice. Stein 
(Ref. 2, p. 288-289] found in his Arctic ice experiments 
that ". . . at no point is the first antisymmetric mode 
[i.e., the flexural wave] supersonic with respect to the 
water; there 1s no coincidence frequency above which 
phase speeds are faster than the water wave speed. This 
implies that in the plane-wave limit the first antisym- 
metric mode never radiates and never loses energy to the 
water." 

After a careful search of recent literature on 
driven plates below coincidence, no theories were found 


that would suggest a radiation pattern other than that 
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of a cosinusoidal pattern. Nor were any theories found 
on the radiation patterns from such acoustically slow 
plates that coincidence was not possible. Therefore, 
theoretical plots, such as those developed for earlier 
sections, were not available. However, some possible 
explanations are offered about the complexities and 
asymmetries of the patterns. 

Since this heed finite-sized plate it could be 
argued that plate modes of vibration were set up, which 
caused nodes in the plate that would affect the pressure 
field in the fluid. However, the pulse mode of opera- 
ihe was used, which reduces the likelihood that this 
was a for discrepancies. Also, it was found 
(and stated earlier in this chapter) that no reflections 
were observed from Ehe Golate sedqe™ a Therefore,s excita- 
tion of normal modes of (iibaticn in the plate are 
thought to be unimportant for the pulse experiment. 

Another possible explanation is that the hydro- 
phone was, at 0.5 m, in the "near-field" of the plate. 
Feit [Ref. 7, p. 1482] used the asymptotic approximation 
koR>>1 in the evaluation of the integral, Equation 3.9. 
At the lowest frequency at wich data were taken, 50 
KHZ, K oR equals 106, so it would appear that the far- 
field assumption is valid. 

The complexities of the patterns might also be 


explained by the pulse not being long enough. The pulse 


-109 


interval was limited to 0.30ms, based on the’ speed of 
the pulse through the water to the receiver and the 
transient response time of the vibrator. Table 3.3 
gives values of travel distance that a flexural wave 
packet would achieve in the plate in 0.30 ms’ for the 


frequencies used. 





The cosinusoidal patterns expected below coinci- 
dence are based on an assumption of cylindrical diver- 
gence of the flexural wave on an infinite plate. That 
is,the pressure at a field point in the fluid is made up 
of the sum of pressures caused by the deflection of each 
incremental area of the plate. If the entire plate were 
not excited, this sum would be incomplete. From Table 
3.2, a pulse of flexural wave of frequency 50 kHz will 
only travel out 0.22m in 0.3 m/s. Thus, the entire 
plate was not being excited. Even at 200 kHz, . the 
entire plate was not excited. This deftainiy could be 
One of the contributing factors in the complexity of 


experimental plots. 
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Prior to the discovery of the electrical pickup 
problem mentioned earlier in this chapter, the plate was 
driven in the CW mode. The complexities in the patterns 
were still there. This would indicate that use of the 
pulse mode was only a contributing factor, not the main 
cause of the complexities of the experimental plots. 

A more important factor is thought to be the ab- 
sorption of the flexural wave in — plexiglass. 
Following the same argument above, that is, the pressure 
at a field point is the sum of the pressures caused by 
the individual deflections of each incremental area of 
the plate, the effect of the flexural wave absorption 
would be to greatly reduce the contributions of radia- 
tions from the incremental areas far from the source. 

The so-called “acoustic mode". of sound propa- 
gation, which this experiment was modelling, assumes 
sound radiation away from the source in the fluid from a 
lossless plate. Plexiglass is a highly absorptive 
medium; therefore, the propagating flexural wave decayed 
in amplitude due to losses as well as divergence as it 
travelled away from the vibrator. Consequently, the 
incremental areas of the plate far from the source did 
not deflect nearly as much as those close to the course, 
so they did not contribute as much to the pressure ata 


given point in the fluid. 
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The plate might be thought of as a complex aper- 
ture, or a quasi-piston, with a radius of the head mass 
radius and a portion of the adjacent plate. A rough 
estimate of the aperture radius is the head mass 
radius plus either one-fourth the flexural wave radius 
plus either one-fourth the flexural wave radius or one 
plate thickness. Both give similar results. Based on 
this, k,a for iil lal into the fluid varies from 2.79 
for a frequency of 50 kHz to 14 at f£ = 250 kHz. 

Even though the aperture is shaded and phased, 
it is noteworthy that for a k,a in excess of 4, a piston 
is beginning to develop a side lobe, and for a koa of 
10, it has about two side lobes. This is not inconsis- 
tent with the angular scale of the lobelike structure 
observed. 

The asymmetry seen in a few of the radiation 
patterns might have been caused by radially non-uniform 
excitation of the plate. As was mentioned in section F 
of this chapter, the vibrator was compressed against the 
plate by a bolt. However, asymmetrical compression of 
the vibrator against the plate was possible. Possible 
causes of asymmetrical compression are improper machin- 
ing of the head mass, imperfect centering of the bolt 


through the vibrator assembly, or a non-uniform layer of 
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grease between the head mass and plate. A non-uniform 
radial excitation of the plate seems to be a plausible 
explanation for asymmetry in the radiation patterns. 

Two important points that can be made are that 
it appears that the radiation pattern is a function of 
the source, not. any flexural wave propagation. 
Secondly, situations where coincidence is not possible, 
such as a plexiglass plate in water, present radiation 
patterns at lower frequencies similar to radiation pat- 
terns below coincidence in situations where coincidence 


is possible. 
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IV. CONCLUSIONS 


This study was a low-cost, "low-tech", first-cut 
attempt to duplicate a sound source mounted on the ice 
in the Arctic. The first experiment performed showed 
that when a plate is introduced at the boundary ofa 
Lloyd's-mirror-type problem, the effect of the plate was 
only to introduce phase delay, causing a shift in the 
node and antinode pattern. If any flexural waves were 
excited in the plate, which might have caused an in- 
crease in amplitude in the received signal, they would 
have been absorbed by the plexiglass. 

In the second experiment performed, it was found 
that even though the theory behind a point-driven fluid- 
loaded plate is well known and many interesting results 
can be gleaned from that theory, extracting ae 
plots upon which a theory can be based is difficult. 

Plexiglass was originally chosen because its speci- 
fic acoustic impedance is fairly close to both ice and 
water. However, the flexural wave absorption suffered 
in commonly available plexiglass and more importantly, 
its low flexural wave speed render it a poor choice to 
use in the lab in this type of experiment, where wave 


generation and propagation are important. 
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Several interesting types of plastics are 
available. Often the manufacturer lists the density and 
Young's modulus of each in its catalog. However, flex- 
ural wave attenuation data are never available. It is 
recommended that the next step in research along these 
lines be to collect samples of each type of plastic and 
measure these parameters. The time Beet lable for this 


study prevented such an analysis. 
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APPENDIX B 
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Figure B.1. Diagram of Tonpilz vibrator 
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